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ABSTRACT 

We investigate a stationary pair production cascade in the outer magnetosphere of 
a spinning neutron star. The charge depletion due to global flows of charged particles, 
causes a large electric field along the magnetic field lines. Migratory electrons and/or 
positrons are accelerated by this field to radiate curvature gamma-rays, some of which 
collide with the X-rays to materialize as pairs in the gap. The replenished charges 
partially screen the electric field, which is self-consistently solved together with the 
distribution functions of particles and gamma-rays. If no current is injected at either 
of the boundaries of the accelerator, the gap is located around the conventional null 
surface, where the local Goldreich- Julian charge density vanishes. However, we first find 
that the gap position shifts outwards (or inwards) when particles are injected at the 
inner (or outer) boundary. Applying the theory to the Crab pulsar, we demonstrate 
that the pulsed TeV flux does not exceed the observational upper limit for moderate 
infrared photon density and that the gap should be located near to or outside of the 
conventional null surface so that the observed spectrum of pulsed GeV fluxes may be 
emitted via a curvature process. Some implications of the existence of a solution for a 
super Goldreich- Julian current are discussed. 



Subject headings: gamma-rays: observation - gamma-rays: theory - magnetic field - 
pulsars: individual (Crab) - X-rays: observation 
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1. Introduction 

The EGRET experiment on the Compton Gamma 
Ray Observatory has detected pulsed signals from 
seven rotation-powered pulsars (e.g., for Crab, Nolan 
et al. 1993, Fierro et al. 1998). The modulation of 
the 7-ray light curves at GeV energies testifies to the 
production of 7-ray radiation in the pulsar magneto- 
spheres either at the polar cap (Harding, Tademaru, 
& Esposito 1978; Daugherty & Harding 1982, 1996; 
Sturncr, Dcrmer, & Michel 1995; Shibata, Miyazaki, 
& Takahara 1998), or at the vacuum gaps in the outer 
magnetosphere (Cheng, Ho, & Ruderman 1986a,b, 
hereafter CHR; Chiang & Romani 1992, 1994; Ro- 
mani and Yadigaroglu 1995; Romani 1996; Zhang 
& Cheng 1997, ZC97). Effective 7-ray production 
in a pulsar magnetosphere may be extended to the 
very high energy (VHE) region above 100 GeV as 
well; however, the predictions of fluxes by the current 
models of 7-ray pulsars are not sufficiently conclu- 
sive (e.g., Cheng 1994). Whether or not the spectra 
of 7-ray pulsars continue up to the VHE region is a 
question that remains one of the interesting issues of 
high-energy astrophysics. 

In the VHE region, positive detections of radiation 
at a high confidence level have been reported from 
the direction of the Crab pulsar (Nel et al. 1993). 
However, as for pulsed TeV radiation, only the upper 
limits have been, as a rule, obtained (Akerlof et al. 
1993; Borione et al. 1997; Srinivasan et al. 1997; 
Yoshikoshi et al. 1997; Sako et al. 2000). If the VHE 
emission originates in the pulsar magnetosphere, a 
significant fraction of it can be expected to show pul- 
sation. Therefore, the lack of pulsed TeV emissions 
provides a severe constraint on the modeling of par- 
ticle acceleration zones in a pulsar magnetosphere. 

In fact, in the CHR picture, the magnetosphere 
should be optically thick for pair-production in or- 
der to reduce the TeV flux to an unobserved level by 
absorption. This in turn requires very high luminosi- 
ties of infrared photons. However, the required IR 
fluxes are generally orders of magnitude larger than 
the observed values (Usov 1994). We are therefore 
motivated by the need to contrive an outer-gap model 
that produces less TeV emission with a moderate in- 
frared luminosity. 

High-energy emission from a pulsar magnetosphere, 
in fact, crucially depends on the acceleration electric 
field, En, along the magnetic field lines. It was Hi- 
rotani and Shibata (1999a, b,c; hereafter Papers I, II, 



III), and Hirotani (2000a; hereafter Paper VI) who 
first considered the spatial distribution of Eu together 
with particle and 7-ray distribution functions. By 
solving these Vlasov equations, they demonstrated 
that a stationary gap is formed around the conven- 
tional null surface at which the local Goldreich- Julian 
charge density, 



Pgj 



2nc 



(1) 



vanishes, where B z is the component of the magnetic 
field along the rotation axis, f2 the angular frequency 
of the neutron star, and c the speed of light. Equa- 
tion (Q) is valid unless the gap is located close to the 
light cylinder, of which distance from the rotation axis 
is given by tzjlc = c/f2. The electrodynamic model 
developed in this paper is basically the same as Pa- 
per VI. However, we find an interesting behavior of 
the gap position, by relaxing the boundary conditions 
to allow electric current injection through the inner 
or the outer boundaries of the gap. 

Subsequently, Hirotani (2000b, hereafter Paper IV; 
2001, Paper V) considered the 'gap closure condi- 
tion' so that a gap may maintain a stationary pair- 
production cascade. In this paper, this closure con- 
dition is generalized into the case when the currents 
are injected through the boundaries. 

In the next two sections, we describe the physical 
processes of pair production cascade and the resultant 
7-ray emission. We then apply the theory to the Crab 
pulsar and present the expected 7-ray spectra in § |[ 
In the final section, we discuss the possibility of a gap 
formation for a super Goldreich- Julian current. 

2. Analytic Examination of the Gap Position 

Let us first consider the gap position analytically 
when there is a current injection into the gap. We 
consider the particle continuity equa tion s in § |2.1| and 
the 7-ray Boltzmann equations in § |2. 



2.1. Particle Continuity Equations 

Under the mono-energetic approximation, we sim- 
ply assume that the electrostatic and the curvature- 
radiation-reaction forces cancel each other in the 
Boltzmann equations of particles. Then the spatial 
number density of the outwardly and inwardly propa- 
gating particles, N+(s) and iV_(s), at distance s from 
the neutron-star surface along the last-open field line, 
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obey the following continuity equations: 
dN ± _ dN± 



dt 



dx 



Q(x), 



where 



1 f°° 

Q(x) = - / de 1 [?7p + G H 
c Jo 



_G_ 



(2) 



(3) 



G+(x, e 7 ) and G_(x, e 7 ) refer to the distribution func- 
tions of outwardly and inwardly propagating 7-ray 
photons, respectively, having energy m e c 2 e 7 . The 
pair production rate for an outwardly propagating (or 
inwardly propagating) 7-ray photon to materialize as 
a pair per unit time is expressed by rj p+ (or r/ p _). 
For charge definiteness, we consider that a positive 
electric field arises in the gap. In this case, N+ (or 
AT_) represents the number densities of positrons (or 
electrons). 

The particle velocity at position (r,6>) becomes 
(eq. [21] in Paper VI) 



(rf2 sin 9 + kB^ 



(4) 



where k is a constant and refer to the azimuthal 
unit vector. In the parentheses, the term rf2sin(9 is 
due to corotation, while kB^ due to magnetic bend- 
ing. Since E\\ arises in the gap, the corresponding 
drift velocity appears as —cE^B^/B 2 . Unless the gap 
is located close to the light cylinder, we can neglect 
the terms containing B^ as a first-order approxima- 
tion. We thus have 



v w v p + rfl sin 6e^. 
Imposing a stationarity condition 

[dt + (rfisin0)fy]iV ± = 0, 



(5) 
(6) 

reminding that the projected velocity on the poloidal 
plane is v p = ccos$i? p /i? p , and utilizing div_B p w 
divi? = 0, we obtain 



±B 



ds\ B 



ccos 



1 f°° 



P+G+ + ?? P _G_ 



(7) 

where $ refers to the projection angle of the particle 
three-dimensional motion onto the poloidal plane. It 
is defined by $ = arcsin(r cnt f2 sin#/c), where r cnt is 
the distance of the gap center from the star center. 
The pair production rate per unit time by a single 
7-ray photon, 77 P ±, are defined as 



v P ±(^) = (1- Mc)c 



dN. 
~ck 



-<T p (e 7 ,e x ,/i c ), (8) 



where er p is the pair-production cross section and 
cos _1 /i c refers to the collision angle between the 7- 
rays and the X-rays (see Paper VI for more details 
about eq. ||); e t h = 2/[(l — // c )e 7 ]. The adopted 
value of /i c will be detailed in § |5.2[ The quantity e x 
refers to the X-ray energy in the unit of m e c 2 . 
Although $ 7^ is adopted after 
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this 

section we simply neglect the projection effect of the 
poloidal velocity and put $ = 0. Then equation (||) 
gives 



± B- 



ds 



N± 
B 



1 r°° 

— / de 7 (G, 
a p Jo 



G_) 
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where G+(s,e 7 ) and G_(s,e 7 ) refer to the distribu- 
tion functions of the outwardly and inwardly propa- 
gating 7-rays; the mean free path A p is defined by 



1 Jo llv+G+de ' ] 



(10) 



Since W <C vulc is justified for the Crab pulsar (Pa- 
per V) , we regard A p to be constant in the gap in this 
section. 

2.2. Boltzmann Equations for Gamma-rays 

Unlike the charged particles, 7-rays do not propa- 
gate along the magnetic field line at each point, be- 
cause they preserve the directional information where 
they were emitted. However, to avoid complications, 
we simply assume that the outwardly (or inwardly) 
propagating 7-rays dilate (or constrict) at the same 
rate with the magnetic field. This assumption gives 
a good estimate when W <C wlc holds. We then 
obtain (Paper VI) 



ds\B 



ccos$ 



G-i 



Vc 



ccos<l> 



where (e.g., Rybicki, Lightman 1979) 



73e 2 r l_ p I 



hRc e-. 



1 3 hcT 3 
m e c 2 47r Rc 



F{s) = s I Ks(t)dt; 



N±, (11) 



(12) 



(13) 
(14) 
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Rc is the curvature radius of the magnetic field lines 
and if 5/3 is the modified Bessel function of 5/3 order. 
The effect of the broad spectrum of curvature 7-rays 
is represented by the factor F(e 7 /e c ) in equation (|l2]). 

Noting that the absorption due to pair production 
is negligible compared with curvature emission term 
on the right-hand side of equation ( pi] ) , and putting 
$ = again, we obtain 



± 



d 



B 



G±(s,e 7 ) 



r?c(e 7 ) N±{a) 
c B 



(15) 



Integrating equation ( |15| ) over e 7 , and combining with 
equation (g), we obtain 



± 



ds 2 



B 



1 N+ — N- 



A p c 



B 



77 c (e 7 )de 7 , (16) 



where /3jv is the upper cutoff dimensionless 7-ray en- 



ergy. In the present paper, we set (3 



N 



= io 5 - 5 



(see§y). 

One combination of the two independent equa- 
tions (|l^) yields the conserved current per magnetic 
flux tube, 



n 
2^ 



Jtot 



ce- 



N+( s)+N_{s) 
B(s) 



(17) 



If itot = 1-0, the conserved current density becomes 
its Goldreich- Julian value. Another combination of 
equations ( |l6| ) gives 



ds 2 



N + - N- 
B 



4N 7 N+ 



B 



where 



W/2 



?7c(e 7 )de 7 



(18) 



(19) 



referes to the expectation value of the number of 7- 
rays emitted by a single particle that runs a typical 
length W/2 in the gap. 

In a stationary gap, the pair production optical 
depth, W/A p , must equal the expectation value for 
a 7-ray to materialize with the gap, A r ~ 1 (iga P /itot)- 
We thus obtain the following condition: 



W 



N 7 jtot ' 



(20) 



which is automatically satisfied by the stationary 
Vlasov equations. Here, the dimensionless current 



density, j gap , created in the gap is defined by 



2nce 



J gap 



_ N + (s 2 ) N+( Sl ) 
B(s 2 ) B(ai) 



B( Sl ) B(a 2 ) 



(21) 



where Si and s 2 designate the position of the in- 
ner and the outer boundaries, respectively. That is, 
W = 82 — si. Equation ( f20| ) corresponds to a general- 
ized verion of the gap closure condition considered in 
Papers IV and V (e.g. eq. [30] in Paper V), in which 
ji = J2 = and hence j gap = jtot was assumed. 
When there is a current injection (i.e., when j\ or j 2 
is non- vanishing) , not only the produced particles in 
the gap but also the injected particles contribute for 
the 7-ray emission. Therefore, the gap width is ad- 
justed smaller compared with ]\ = j 2 = case by 
the factor j gap /jtot- Utilizing condition (pfj|), we can 
rewrite equation Jig) into the form 



ds 2 



B 



= A 2^± N + - N -, (22) 

JgapTU B { ' 



To solve the differential equation (P2|), we impose 
the following two bounday conditions: 



C6 B( Sl ) 2^ Jl ' 

JV-(a 2 ) 0. . 
ce „, \ = — j 2 - 



(23) 



(24) 



B(a 2 ) 2vr- 

With the aid of equation (pi]), these two bounday 
conditions give 

5 = -^-Ugap ~ h+32) (25) 

at s = Si, and 



= ^-C/gap +h - h) (26) 



/; 2w u 

at s = S2- It follows from equation (|2 



N + - N. 
B 



2nce 



sinh 



Jtot a - S cnt 
Jgap W/2 



Sim, ( 4 /ST) 

J gap / 



cosh 



+ (ji -ia)- 



Jtot a — S cn t 
jgap W/2 



cosh 



Jtot \ 
Jgap / 



(27) 
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where the gap center position is defined by 

Sent — 2 

2.3. Poisson Equation 



(28) 



The real charge density e(N + — AT_), which is given 
by equation (|2^) appears in the Poisson equation for 
the non-corotational potential \P. Neglecting rela- 
tivistic effects, and assuming that typical transfield 
thickness of the gap, D±, is greater than or compara- 
ble with W, we can reduce the Poisson equation into 
the one-dimensional form (Paper VI; see also § 2 in 
Michel 1974) 



V 2 * = 4tt 



e(N+ - N-) + 



nB 7 



2ttc 



(29) 



where e designates the magnitude of the charge on 
an electron. Substituting equation (|27| ) into (p9|), we 
obtain 



V 2 * 



2BQ 



Jgap/odd 



+ (jl -32) fe- 



W/2 
S - Sent \ Bz 

W/2 ) B 



where 



sinh 



/odd (a) 



Jtot 
3 gap 



■ X 



and 



sinh 



cosh 



/even(^) 



Jtot 
Jgap 

jtot 
Jgap 



, (30) 



(31) 



cosh 



Jtot 
J gap 



(32) 



There are essentially three assumptions that are used 
to derive equations (30), ( |3l] ) , and ( p2|) : the radiation- 
reaction forces exactly cancel with the electrostatic 
force in the particles' Boltzmann equations; ?7 P +(e 7 ) = 
?7 P _ (&y), which may be justified for a power-law, mag- 
netospheric X-ray component; and the Poisson equa- 
tion is analyzed one-dimcnsionally along the magnetic 
field line. 

2.4. Generalization of the Null Surface 

To examine the Poisson equation ( |30| ) analytically, 
we assume that the transfield thickness of the gap 
is greater than W and replace V 2 * with cP^f/ds 2 . 



Furthermore, we neglect the current created in the 
gap and simply set j gap = 0. 

First, consider the case when a current injects from 
neither of the boundaries, that is, j\ = j% = 0. It fol- 
lows that the derivative of the acceleration field (i.e., 
—d 2 ty/ds 2 ) vanishes at the conventional null surface 
where B z vanishes. We may notice that —d 2 ^/ds 2 is 
positive at the inner part of the gap and changes its 
sign near the gap center (s = s cn t) to become nega- 
tive at the outer part of the gap. Therefore, we can 
conclude that the gap is located (or centers) around 
the conventional null surface, if there is no current 
injection from outside. 

Secondly, consider the case when a current is in- 
jected at the inner boundary (at s = Si) and j\ —ji > 
holds. Since the function / eve n is positive at arbi- 
trary s, the gap center is located at a place where 
B z is negative, that is, outside of the conventional 
null surface. In particular, when j\ — j ' 2 ~ 1 holds, 
—d 2 "^ /ds 2 vanishes at the place where B z ~ — B. In a 
vacuum, static dipole field, B z ~ — B is realized along 
the last-open field line at the light cylinder. There- 
fore, the gap is expected to shift towards the light 
cylinder, if the injected current density at the inner 
boundary approaches the Goldreich- Julian value. We 
may notice here that / even is less than unity, because 
\s — s cn t I does not exceed W/2. 

Thirdly and finally, consider the case when j\ — 
j 2 —l holds. In this case, —d 2x S/ds 2 vanishes at the 
place where B z ~ B. Therefore, gap is expected to be 
located close to the star surface, if a Goldreich- Julian 
current density is injected at the outer boundary. In 
what follows, we will examine more accurately these 
predictions on the gap position vs. current injection, 
by solving the Vlasov equations (0), (|ll|), and (30) 
numerically. 

3. Basic Equations and Boundary Conditions 

In the present paper, we assume that the transfield 
thickness, D±, of the gap is much greater than W, 
and neglect the transfield derivatives in the Poisson 
equation ( |29| ) . We consider that this one-dimensional 
analysis could be justified because D± ~ 6W is re- 
quired so that the predicted GeV flux may be con- 
sistent with the EGRET observations (§ |5.3.l|) . We 
rewrite the Vlasov equations into the suitable forms 
for numerical analysis in § 3.1, and impose boundary 



conditions in S 3.2 
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3.1. One-dimensional Vlasov Equations 

As will be shown at the end of this section, it is 
convenient to introduce the typical Debey scale length 



' W riB cnt 

m c 2nce ' 



(33) 



where B cnt represents the magnetic field strength at 
the gap center. The dimensionless coordinate variable 
then becomes 

C = (wp/c)a. (34) 

By using such dimensionless quantities, we can rewrite 
the Poisson equation into 



KK)-n-K)] + 



(35) 



(36) 



d£, B cn t B cn t 

where = e^(s)/ (m e c 2 ); the particle densities per 
unit flux tube are defined by 



n±(0 



2-Kce N ± 



(37) 



We evaluate B z /B at each point along the last-open 
field line, by using the Newtonian dipole field. 

Let us introduce the following dimensionless 7-ray 
densities in the dimensionless energy interval between 
ft_i and ft: 

1-nce l' l3i 

9±(0=7^- rfe T G±( S ,e 7 ). (38) 

"-Dent J/3i_i 

In this paper, we set ft — 10 2 , which corresponds 
to the lowest 7-ray energy, 51.1 MeV. We divide the 
7-ray spectra into 9 energy bins and put ft = 10 25 , 
ft = 10 3 , ft = 10 3 - 5 , ft = 10 4 , ft = 10 4 - 5 , ft = 
10 4 - 75 , ft = 10 5 . ft = 10 5 - 25 , and ft = 10 5 - 5 . 

We can now rewrite the continuity quation (0) of 
particles into 

i=l 

where the magnetic field strength, B, is evaluated at 
each £. The dimensionless redistribution functions 
7? p± are evaluated at the central energy in each bin as 



1 

V P ± = —v P ± 



ft-i + ft 



(40) 



A combination of equations ( J39| ) gives the current 
conservation law, 



j tot = n+(£) + n_(f) = constant for £, 



(41) 



which is equivalent with equation (gj 

The Boltzmann equations ( |ll| ) for the 7-rays are 
integrated over e 7 between dimensionless energies 
ft_i and ft to become 



where i = 

4 



d£ v ' cos$ * B cnt cos $ 



n±, (42) 



1, 2, • • •, m (to = 9) and 
V3e 2 r 



Lu p hR c 



ds 



Kb (t)dt (43) 



is dimensionless. 

Equating the electric force e\d*k /dx\ and the radi- 
ation reaction force, we obtain the saturated Lorentz 
factor at each point as 



Isat — 



/3i? c 



V 2e 



ds 



1 



1/4 



(44) 



we compute the curvature radius R c at each point for 
a Newtonian dipole magnetic field. Since the maxi- 
mum of \d*5f /dx\ and the potential drop are roughly 
proportional to W 2 and W 3 , respectively (Paper V), 
the particles become unsaturated for very small W. 
To avoid an overestimation of the Lorentz factor in 
such cases, we compute L by 



1 



1 



1 



r Vr-t 2 ^ 2 (6)' 



(45) 



where ^(£2) represents the maximum attainable Lorentz 
factor. 

3.2. Boundary Conditions 

We now consider the boundary conditions to solve 
the Vlasov equations ©, @, ©, and ©. At 
the inner (starward) boundary (£ = £1), we impose 
(Paper VI) 

£||(6) = 0, (46) 

= 0, (47) 

4(6) =0 (» = 1,2,...,9). (48) 

It is noteworthy that condition ( pfSj ) is consistent with 
the stability condition at the plasma- vacuum interface 
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if the electrically supported magnetospheric plasma is 
completely-charge-separated, i.e., if the plasma cloud 
at £ < £1 is composed of electrons alone (Krause- 
Polstorff & Michel 1985a,b; Michel 1991). We as- 
sume that the Goldreich-Julian plasma gap boundary 
is stable with E\\ — on the boundary, £ = £1. 

Since positrons may flow into the gap at £ = £1 as 
a part of the global current pattern in the magneto- 
sphere, we denote the positronic current per unit flux 
tube at £ = £1 as 



which yields (eq. [ftl"|) 



= jtot ~ji 



(49) 



(50) 



At the outer boundary (£ = £ 2 )j we impose 

^||(6) = 0, (51) 



5_(6) = (t = l J 2,..-,9) J 



(52) 
(53) 

Conditions fl49| ) and (p^) are equivalent with (|2^) and 
©• 

The current density created in the gap per unit flux 
tube can be expressed as 



Jgap — Jtot — ]i — 32- 



(54) 



This equation is, of course, consistent with equa- 
tion (^TJ). We adopt j gap , j\, and j2 as the free pa- 
rameters. 

We have totally 24 boundary conditions (|i^)~(|5^) 
for 22 unknown functions En, n±, g± (i — 1,2, • • 
•,9). Thus two extra boundary conditions must be 
compensated by making the positions of the bound- 
aries £i and £2 be free. The two free boundaries ap- 
pear because E\\ = is imposed at both the bound- 
aries and because j gap is externally imposed. In other 
words, the gap boundaries (£1 and £2) shift, if ji 
and/or j% varies. 

Let us briefly comment the convenience of the in- 
troduction of the dimensionless coordinates and vari- 
ables. It follows from the Vlasov equations (|35|), d36|), 
(39), and ( [i"2| ) that the solutions £, ip, Eu , n±, and g l ± 
are unchanged if B, ry p± , and rf c are invariant. Con- 
sider the case when the normalization of dN x /de x is 
doubled. In this case, equations (Sh and (|4(]) show 
that rjp is invariant if we also double w p . Note that e c 



in equation (43) is proportional to w p ^ 4 rV 2 (eq. fl3]| ), 
where T - r sat cx Rl /2 \d^ /ds] 1 / 4 is used. It follows 
that e c is also invariant if we increase R c by 2 -3 / 2 
times. It should be noted that 77* (eq. E|) is invari- 
ant by this change of parameters. On these grounds, 
we can reduce one degree of freedom in the free pa- 
rameters. 

4. Predicted Gamma-ray Flux 

In this section, we detail the method how to com- 



pute vF v spectrum in GeV energies in 
TeV energies in 



4.2 



4.1 



and in 



4.1. GeV Spectra 

The GeV spectra of outwardly and inwardly prop- 
agating 7-rays are obtained from #+(£2) and 
At position £, the 7-ray emission rate becomes 

O Fi 

7 -ray flux = A cr (£) c ~^9±(0 (55) 

where A CI refers to the cross section of the gap at £. 
Multiplying the mean 7-ray energy */ /3i/?i_im e c 2 , on 
equation J55|), dividing it by AflQ C yd 2 , and further 
dividing by the frequency interval m e c 2 (/3i — P\-x)/h, 
we obtain the flux density, F v ; here, A^GcV is the 
emission solid angle, and h the Planck constant. We 
thus obtain the GeV flux 



vF v = 



AA-l 2 

-m c c 



riB r 



■fj± 



A - &- 



2ire Ail 



GeV 



d 2 ' 



(56) 



To compute the 7-ray flux emitted outwardly (or in- 
wardly) from the gap, we adopt the plus (or the mi- 
nus) sign in g± and ev alua te A cr g + at £ = £2 (or £1). 
As will be shown in § 5.3, W <C ^lc holds for the 



Crab pulsar. We thus simply apply the same cross 
section for both the outwardly and inwardly emit- 
ted 7-rays and put A CI = D\, where D± should be 
greater than or at least comparable with W for the 
one-dimensional approximation of the Poisson equa- 
tion (§ |2^) to be justified. 

It is noteworthy that the particles lose most of their 
energy in the gap if / acc *C W holds, where l acc refers 
to the length scale for particles to be accelerated to 
the saturated Lorentz factor (eq. 0]). That is, we 
can neglect the primary luminosity emitted by the 
particles running outside of the gap, compared with 
that emitted by the particles running inside of the 
gap, if ^ acc <C W. Since the mono-energetic approxi- 
mation of the particle motion (§ |2.l|) is justified when 
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^acc "C W, the neglect of GeV emission by the par- 
ticles running outside of the gap is consistent with 
the mono-enegetic approximation. We thus compute 
the GeV luminosity from the solved 7-ray distribution 
functions g\{£,2) and <7!_(£i). 

4.2. TeV Spectra 

Once the electrodynamic structure of the gap is 
solved, we can further compute the upscattered 7-ray 
flux emitted from the whole accelerator, if addition- 
ally give the infrared photon field. This treatment is 
justified unless the upscattered, TeV luminosity ex- 
ceeds the curvature-radiated, GeV one. 

If an electron or a positron is migrating with 
Lorentz factor r > 1 in an isotropic photon field, 
it upscatters the soft photons to produce the follow- 
ing number spectrum of 7-rays (Blumenthal & Gould 
1970): 



dN 
dtde^ 



3 c dN m de m 
-<7t; 



4 T 2 de m ei R 



2qlnq+{l + 2q)(l-q) + 



(Qq) 2 (l-q) 



2(1 + Qq) 



,(57) 



where Q = 4erR.r, q = e 7 /Q(T — e 7 ), dNm/dem refers 
to the IR photon density per unit dimensionless enegy 
interval between eiR and eiR + dem, ctt is the Thom- 
son cross section; eiR and e 7 are the energies of the IR 
and the upscattered photons in units of m c c 2 . Equa- 
tion ( |57j ) is valid if the resonance effects are negligible, 
that is, B < B crit = 4.4 x 10 13 G. This inequality is 
satisfied except for the polar cap. The flux density of 
the upscattered photons becomes 



TP - N ° h / ' £IR ' 

v ~ Arixcv d 2 ' 7 



dN 



dtde-ydtiR, 



dem. (58) 



where AfixeV refers to the emission solid angle of the 
upscattered photons. In this paper, we estimate N e 
with 



Ac = (jgap 



2irce 



-WD 1 



(59) 



to compute the outwardly propagating TeV flux, 
which are emitted by outwardly propagating particles 
(i.e., positrons) and with 



(60) 



to compute the inwardly propagating TeV flux, which 
are emitted by inwardly propagating particles (i.e., 
electrons). 



Multiplying the 7-ray frequency e 7 m c 2 //i on the 
F v flux density (eq. |Q), we obtain the upscattered 
flux 



(jgap + Ja)m e c : 



flB r 



2ne An TcV d 2 7 
\ir / ... . dent, (61) 



em, 1. 



cdtde^deiYi 



where j a = j± (or j'2) for outwardly (or inwardly) 
emitted 7-rays. As the emission solid angles, we as- 
sume 

2nW 

Ail GcV = AOtoV = (62) 

in this paper. 

We also consider the extrinsic absorption of the 
TeV photons outside of the gap. For a homogeneous 
and isotropic IR field, the optical depth becomes 



r(e 7 ) 



dN 



IR 



dc 



IB 



fp(eiR, e 7 , ^ c )rfeiR, 



(63) 

where the path length is assumed to be tn LC /2. We 
apply the same path length for all the cases considered 
so that the extrinsic absorption may work equally. 

5. Application to the Crab Pulsar 
5.1. Input Infrared Field 

Consider the case when the IR spectrum is homo- 
geneous and expressed by a single power-law, 



dNm 



= N e% 



IR, 



(64) 



where No and a are spatially constant. For an 
isotropic field, the specific intensity becomes 



I v = —hN^ a+1 



An 



^ir 



(65) 



Assuming that this uniform sphere has radius wlc, 
we obtain the following flux density at distance d: 

K . f(^2)W+- ^ 

= 4.5 x m-^n,- 2 (jjj) Woem" +1 Ji66) 

As the lower and upper cutoff IR photon energies, 
we adopt eiR, m i n = 10" 8 and eiR, ma x = 10~ 6 , where 

ElR.min < ^IR < ClR.max- 
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Because the pulsed flux around eV energies are dif- 
ficult to be observed, we consider the following two 
cases for the set of Nq and a: 

case A We assume that the IR spectrum below 
em < 1CT 6 (or equivalently, below 1.23 x 10 14 Hz) 
are optically thick for synchrotron self-absorption and 
adopt a — 1.5. Setting F v — 3 mJy at eiR = 10~ 6 , 
which is consistent with near-IR and optical obser- 
vations (Eikenberry et al. 1997), we obtain No = 
2.3 x 10 32 cm- 3 . 

case B Interpolating the phase-averaged color spec- 
trum in UV, U, B, V, R (Percival et al. 1993), J, H, 
K (Eikenberry et al. 1997) bands, and the radio ob- 
servation at 8.4 GHz (Moffctt and Hankins 1996), we 
obtain N = 1.5 x 10 17 cm" 3 and a = -0.88. In fig- 
ure ||, we present the fitted spectrum with the solid 
line; the ordinate is vF v in Jy Hz. 

5.2. Input X-ray Field 

HEAO 1 observations revealed that the X-ray spec- 
trum in the primary pulse phase is expressed by 



dN pl 
de x 



N p ie x a (e n 



< e x < en 



(67) 



witha = -1.81andiV p i = 5.3xl0 15 (d/kpc) 2 (r cnt /cc7 LC )~ 2 
(Knight 1982). We adopt e min = 0.1keV/511keV and 
e max = 50keV/511keV. Unlike the IR field, which is 
assumed to be homogeneous within radius va^c, we 
suppose that the X-rays are emitted near to the gap. 
In this case, the X-ray density computed from the ob- 
served flux will increase as the gap is located close to 
the star. To consider such effects, we simply assume 
that the X-ray density is proportional to the inverse 
square of r cnt . 

The angle dependence of the specific intentity of 
the X-ray field is considered in the collision angle, /i c 
(eq. ||, or ||t|). In the case of the Crab pulsar, the 
X-ray field is dominated by a power-law component, 
which is probably emitted near the outer-gap accel- 
erator rather than from the neutron star surface. We 
thus simply evaluate the cosine of the collision angles 
as 

fi c = cos{W/zu LC ), (68) 

for both inwardly and outwardly propagating 7-rays. 
Aberration of light is not important for this com- 
ponent, because both the X-rays and the 7-rays are 
emitted nearly at the same place. We may notice here 
that this is a rough estimate of /j, c and that e t h = 
2/[(l — ^0)^7] strongly depends on ^ c if W <C w\,c 



(i.e., if 1 — ^ c <C 1). In the case of the Crab pulsar, 
W/zuhc ~ 0.05 holds (see fig. 2); therefore, the true 
results will depend on the detailed beaming geometry 
of the secondary X-rays, which are emitted outside 
of the gap along local magnetic field lines via syn- 
chrotron process. However, to inquire into this mat- 
ter would lead us to into that specialized area of the 
magnetic field configuration close to the light cylin- 
der. Such a digression would undoubtedly obscure the 
outline of our argument. 

5.3. Results 

Let us now substitute the X-ray field into equa- 
tion (||) and solve the Vlasov equations by the method 
described in § 2. It should be noted that we do not use 

to solve £j|(£): 



the IR spectrum considered in § 5.1 



n ±(£,), an d <7± z (£); they are necessary when we con- 
sider the TeV emission due to IC scatterings. 

In Paper I, it is argued that the IC scatterings dom- 
inates the curvature radiation in the outer gap in the 
Crab pulsar magnetosphere if we adopt /i c — 0. How- 
ever, in the present paper, we consider that the col- 
lision angles are much smaller than 90° and adopt 
/i c = cos(M^/n7Lc)- In this case, unless the gap is lo- 
cated well inside of the conventional null surface, cur- 
vature radiation becomes the primary process in 7-ray 
emission, and hence in the radiation-reaction forces. 
The rotational frequency and the magnetic moment 
are Q = 188.1rad s" 1 and n m = 3.38 x 10 30 G cm 3 . 

5.3.1. Electric Field Structure 

To reveal the spatial distribution of the accelera- 
tion field, we consider four representative boundary 
conditions: 

case 1 {31,32) = (0,0) — * solid curves, 
case 2 {31,32) = (0.3,0) — * dashed curves, 
case 3 {31,32) = (0.6,0) — * dash-dotted curves, 
case 4 {31,32) = (0,0.3) — + dotted curves. 
That is, for case 2 (or case 4), the positronic (or elec- 
tronic) current density flowing into the gap per unit 
flux tube at the inner (or outer) boundary is 30% 
of the typical Goldreich- Julian value, il/2n. We fix 
jgap = 0.01 for all the four cases, because the solu- 
tion forms a 'brim' to disappear (fig. 2 in Hirotani 
& Okamoto 1998) if j gap exceeds a few percent. In 
what follows, we adopt 45° as the magnetic inclina- 
tion, which is necessary to compute B at each point 
for the Newtonian dipole field. Since the X-ray field is 
dominated by a power-law component, the inclination 



9 



does not affect /x c . 

The resufts oiE\\ (£) for the four cases are presented 
in figure |^. The abscissa designates the distance along 
the last-open field line and covers the range from the 
neutron star surface (s = 0) to the position where the 
disance equals s — 1.2tz7lc = 1-91 x 10 6 m. 

The solid curve (case 1) shows that the gap is lo- 
cated around the conventional null surface. However, 
the gap shifts outwards as j\ increases, as the dashed 
(case 2) and dash-dotted (case 3) curves indicate. 
This result is consistent with what was predicted in 
Shibata and Hirotani (2000) analytically. 

On the other hand, when j 2 increases, the gap 
shifts inwards and the potential drop, \I/(s2), re- 
duces significantly. For example, we obtain ^(s2) = 

7.1 x 10 12 V for case 4, whereas 1.7 x 10 13 V for case 2. 
A detailed physical interpretation is given in § |6.1| . 

5.3.2. Gamma-ray Spectra 

We compute the GeV and TeV spectrum by the 
method described in § |]. We adopt the cross sec- 
tional area of D± 2 = (6W) 2 for all the cases to be 
considered, so that the GeV flux in cases 1 and 2 may 
be consistent with observations. If D±_ increase twice, 
both the GeV and TeV fluxes increases four times. 

First, we consider case A in which the IR spectrum 
is approximated by SSA with turnover frequency ~ 

1.2 x 10 14 Hz. In this case, the pair-production opti- 
cal depth r computed from equation ( |63|) becomes as 
presented in figure [| This result indicates that the 
TeV flux is significantly absorbed above 1 TeV. 

For the four different boundary conditions (cases 1, 
2, 3, and 4), we present the spectra of the outwardly 
and inwardly propagating 7-rays in figures ^ and |^, re- 
spectively. In GeV energies, the observational pulsed 
spectrum is obtained by EGRET observations (open 
circles; Nolan et al. 1993), while in TeV energies, 
only the upper limits are obtained by Whipple obser- 
vations (open squares; Weekes et al. 1989; Reynolds 
et al. 1993; Goret et al. 1993; Hillas, A. M.; Lessard 
et al. 2000), Durham observations (open triangle; 
Dowthwaite et al. 1984), and CELESTE observations 
(open square at 60 GeV; Holder, J., private commu- 
nication). The filled circles denote the unpulsed flux 
obtained by CANGAROO observations (Tanimori et 
al. 1998). 

It follows from figures || and || that the TeV flux is 
undetectable except for hv ~ 10 TeV. Around 10 TeV, 
the 7-ray flux is slightly less than or comparable with 




frequency [Hz] 



Fig. 1. — A single power-law fit of phase-averaged 
color spectrum of the Crab pulsar (case B, see § 5.1). 
The abscissa is the photon frequency in Hz, while the 
ordinate is the photon flux in Jy Hz. 



Crab 



inclination= 45 deg 




5x10 J 10° 

distance along field line [meter] 

Fig. 2. — Distribution of E\\(s) for the Crab pulsar 
with ot\ = 45°; the abscissa is in meters. The solid, 
dashed, dash-dotted, and dotted curves correspond to 
the cases 1, 2, 3, and 4, respectively (see text). 
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the observational upper limits for cases 1, 2, and 3, 
and exceeds the limits for case 4. Nevertheless, we 
can exclude case 4 from consideration, because the 
expected GeV spectrum is very very soft and is in- 
consistent with the EGRET observations, whatever 
D±_ we may assume. 

It is noteworthy that the GeV spectrum, which 
does not depend on the assumed IR field, depends 
on ji and ji significantly. In particular, in case 4 (as 
the dotted curves show), the GeV emission signifi- 
cantly decreases and softens, because both the poten- 
tial drop and the maximum of Bn reduce as the gap 
shifts inwards. As a result, it becomes impossible to 
explain the EGRET flux around 10 GeV, if the gap 
is located well inside of the conventional null surface. 

Next, let us next consider the case B in which the 
IR spectrum is interpolated from radio and optical 
pulsed fluxes. In this case, the pair-production op- 
tical depth r computed from equation( ^33|) becomes 
as presented in figure [| Therefore, the emitted TeV 
flux significantly reduces above 1 TeV. 

The spectra of the outwardly and inwardly prop- 
agating 7-rays are presented in figures and ||, re- 
spectively. It follows from the two figures that the 
TeV fluxes exceed the observational upper limits for 
cases 2, 3, and 4. In case 1, the upscattered flux is 
small because of its small N c , which is proportional 
to j tot = 0.01. 

Let us briefly consider the case when the interpo- 
lated spectrum (No = 1.5 x 10 17 and a = —0.88) ex- 
tends to much higher frequencies and adopt eiR. max = 
10" 5 (or 1.2 x 10 15 Hz). In this case, r > 2 holds 
above 0.2 TeV (fig. §); the absorbed TeV flux is thus 
suppressed below the observational upper limits as 
figure [l(] indicates. 

In short, we can conclude that the problem of the 
excessive TeV flux does not arise if the IR field is 
represented by a SSA spectrum (case A) or if the IR 
field is interpolated by a single power law (case B) 
with a large cut-off energy (eiR imax ~ 10~ 5 ). 

5.4. Dependence on Magnetic Inclination 

To investigate how the results depend on the mag- 
netic inclination, we present the expected Crab pul- 
sar spectra for ot\ — 75° in figure [n|. The dashed 
and dash-dotted curves correspond to cases 2 and 
3, while the dash-dot-dot-dotted ones to the case of 
ji = 0.5 and j 2 = 0. Case 1 is not depicted because 
the central energy of curvature-radiated photons be- 



Crab 



co - 




photon energy [eV] 



Fig. 3. — Pair-production optical depth for the IR 
field represented by SSA spectrum (case A). 



Crab 

inclination= 45 deg 




photon enerqy TeV] 

Fig. 4. — Spectra of the outwardly propagating 7- 
rays emitted from the Crab pulsar magnetosphere 
when the IR field is approximated by a SSA spectrum 
(case A). The solid, dashed, dash-dotted, and dotted 
curves correspond to the same boundary conditions 
as in figure 2. 
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Crab 

inclination= 45 deg 




photon enerqy [eV] 

Fig. 5. — Same figure as figure 4 but the 7-rays are 
inwardly propagating. 



Crab 



co - 

CD CO - 




photon energy [eV] 



Fig. 6. — Pair-production optical depth for the single 
power-law IR spectrum (fig. 1) fitted between radio 
and optical bands (case B). 



Crab 



inclination= 45 deg 




photon enerqy [eV] 

Fig. 7. — Spectra of outwardly propagating 7-rays 
from the Crab pulsar magnctosphere when the IR 
spectrum is interpolated from radio and optical bands 
with a single power law (case B). The solid, dashed, 
dash-dotted, and dotted curves correspond to the 
same cases as in figure 2. 
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inclination= 45 deg 




photon enerqy [eV] 

Fig. 8. — Same figure as figure 6 but the 7-rays are 
inwardly propagating. 
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Crab 




10 



10" 10' 
photon energy [eV] 



Fig. 9. — Same figure as fig. 6 but with a large upper 
cutoff energy, eiR, max = 1(T 5 . 



Crab 




nclination= 45 deg 



D_perp= 6 Wi 



OL 



photon enerqy [eV] 

Fig. 10. — Gamma-ray spectra from the Crab pulsar 
magnetosphere for case B but with a large upper cut- 
off energy, eiR jmax = 1CP 5 . The thick (or thin) curves 
denote outwardly (or inwardly) propagating 7-rays. 



comes comparable with /?nm c c 2 = 90.8 GeV; in this 
case, its hard spectrum would be inconsistent with 
the EGRET pulsed spectrum below 30 GeV and the 
CELESTE upper limit at 60 GeV. Moreover, case 4 
(i.e., ji = and ji = 0.3) is excluded in figure |TT]; this 
is because the gap is located so close to the star sur- 
face that the IC scatterings dominates the curvature 
process. 

Comparing the GeV spectra in figure [ll] with those 
in figures ^ and ^, we can confirm that the curvature 
emission becomes hard and luminous if a\ increases. 
Its physical interpretation will be discussed in § |6.2. 



In figure [ll], case A is adopted as the infrared spec- 
trum; however, the IR field is not important when we 
discuss the curvature-radiated 7-ray spectrum. 

It also follows from figure [ll] that the observed, 
pulsed GeV spectrum can be explained if we take 
ji = 0.5 and j% = for a; = 75°. In other words, 
the curvature spectrum becomes analogous between 
ji = 0.3 for cti = 45° (dashed curve in fig. [|) and 
ji = 0.5 for «i = 75° (thick, dash-dot-dot-dotted one 
in fig. [ll]), if we fix j gap = 0.01 and j'2 = 0. That 
is, a greater j\ is preferable for a greater a{. It is 
natural, because the decrease of the distance of the 
intersection between the conventional null surface and 
the last-open field line from the star surface (with in- 
creasing ai) should be compensated by shifting the 
gap outwards (with increasing ji), so that the gap 
may have the comparable magnetic and X-ray field 
strengths. 

On these grouds, we can conclude that we can- 
not decouple the effects of the magnetospheric cur- 
rents (j'ij'2) and cti, if we only compare the fluxes of 
the outwardly propagating 7-rays (when ji > j'2). It 
would be possible to argue that these two effects could 
be decoupled if we considered the inward/outward 
flux ratio, or the three-dimensional structure of the 
accelerator. However, such details are irrelevant to 
the main subject of this paper. 

6. Discussion 

In summary, we have developed a one-dimensional 
model for an outer-gap accelerator in the magneto- 
sphere of a rotation-powered pulsar. When a mag- 
netospheric current flows into the gap from the outer 
(or inner) boundary, the gap shifts inwards (or out- 
wards). In particular, when a good fraction of the 
Goldrcich- Julian current density is injected from the 
outer boundary, the gap is located well inside of the 
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conventional null surface; the resultant GeV emission 
becomes very soft and weak. Applying this method 
to the Crab pulsar, we find that the gap should be 
located near to or outside of the conventional null 
surface, so that the observed GeV spectrum of pulsed 
GeV fluxes may be emitted via a curvature process. 
By virtue of the absorption by the dense IR field in the 
magnetosphere, the problem of excessive TeV emis- 
sion does not arise. 

6.1. Gap Width vs. Current Injection 

By utilizing the gap closure condition (ptj|), we can 
interpret why W becomes significantly small when the 
gap is located well inside of the conventional null sur- 
face (§ 5.3.1). First, the X-ray density becomes large 
at small radii to reduce A p in equation (p0|). Sec- 
ondly, the ratio jgap/jtot decreases as fa increases. 
As a result, W decreases very rapidly with increasing 
fa. When W decreases, A" 7 decreases to some extent; 
however, this effect is passive and cannot change the 
conclusion. On these grounds, the gap width signif- 
icantly decreases when particles are injected at the 
outer boundary. Therefore, the potential drop also 
decreases significantly. 

On the other hand, when the gap is located out- 
side of the conventional null surface, the decreased 
Jgap/jtot due to the increase of fa partially cancels 
with the increase of A p due to the diluted X-ray field. 
Thus, the gap width is roughly unchanged when par- 
ticles are injected at the inner boundary. 

6.2. Interpretation of the Magnetic-Inclination 
Dependence 

In this subsection, we interpret the dependence of 
the results on a; (in § 5^). In Pape V, it was predicted 
that W (= 2H in their notation) is a decreasing func- 
tion of «i for all the twelve pulsars considered. The 
reasons are five folds: 

• With no current injection (i.e., fa = fa = Q as con- 
sidered in Paper V), the gap is located at the inter- 
section of the last-open field line and the conventional 
null surface, where B z vanishes. 

• The intersection approaches the star if a; increases. 

• The density of the X-ray field illuminating the gap 
increases as the intersection approaches the star (or 
equivalently, as r cnt decreases). 

• It follows from the closure condition (eq. [^0|) that 
W ex A P /7V 7 . If we neglect the variations in N 7 , W is 



proportional to A p cx N K 



7 ' 

1 cx r cnt ~ 2 , where N x is the 



X-ray number density. Therefore, W decreases with 
decreasing r cnt , and hence with increasing a{. 

• In reality, A 7 decreases if W decreases. As a result 
of this 'negative feedback effect', the decrease of W 
for an increasing a\ is partially canceled. However, 
this effect is passive; therefore, the conclusion of the 
decreasing W with increasing a; is unchanged. 

In Paper V, the GeV emission is predicted to be- 
come hard and luminous, as ct\ increases for the same 
set of jgap, ji, and fa. The reasons are fivefold: 

• The gap approaches the star (i.e., r cnt decreases), 
as cti increases for fixed fa and fa (say, fa = fa = as 
considered in Paper V). 

• The magnetic field in the gap increases as B cx 
?~cnt -3 when the gap approaches the star. 

• As a result of this rapid increase of B, E\\ increases 
(e.g., eq. [§9|), in spite of the decreasing W, as stated 
in the paragraph just above. 

• The increased E\\ for a larger ai results in a harder 
curvature spectrum in GeV energies. 

• The potential drop in the gap is roughly propor- 
tional to the maximum of E\\ in the gap times W . 
Because of the 'negative feedback effect' due to A 7 , 
the weakly decreasing W cannot cancel the increase 
of E». As a result, the potential drop, and hence the 
GeV luminosity increases with increasing aj. 

6.3. Super Goldreich-Julian Current 

In this subsection, we briefly discuss the relaxation 
of the limit of the current density flowing in the gap 
along the field lines. In Papers I, II, III, VI, in which 
ji = fa = is assumed, stationary gap solutions were 
found only for a small j ga p- By the revised method 
presented in this paper, the solutions for the Crab pul- 
sar exist for j gap < 0.0255, if we set fa — fa = 0. The 
solution of -E||(s) for fa = fa = and j gap = 0.0255 
for the Crab pulsar when ct\ — 45° is depicted in fig- 
ure [l|. Because of the 'brim' at the inner boundary, 
no solution exists for j gap > 0.0255. In this case, 
itot = jgap + ji + fa is limited only below 0.0255, 
which is much less than the typical Goldreich-Julian 
value, 1. 

Let us briefly consider how much j t ot is needed for 
the observed spin-down luminosity to be emitted. If 
we assume that all the current flowing in the mag- 
netosphere penetrate the gap, then the net current 
becomes J = (il/2ir)j to t'&, where ^ is the magnetic 
fluxes along which the current is flowing. Assuming 
a magnetic dipole radiation, we obtain the potential 
drop at the stellar surface as V* ~ Vl^ /(ttc). The 
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spin-down luminosity then becomes 
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Fig. 11. — Gamma-ray spectra from the Crab pul- 
sar magnetosphere with cti = 75°, when the IR field 
is approximated by a SSA spectrum (case A). The 
thick (or thin) curves denote outwardly (or inwardly) 
propagating 7-rays. 
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Fig. 12. — Distribution of E^(s) for the Crab pulsar 
with «j = 45° when j gap = 0.0255, j\ = 32 = 0. 
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J x K 



Jtot 
C 



(69) 



If the gap is geometrically thick in the transfield di- 
rections, we may expect that the field lines thread 
the polar cap with area A po \ c = 7r(r* sin #*) 2 , where 
r* refers to the stellar radius and 6*„ to the colatitude 
angle between the magnetic axis and the last-open 
field line. Utilizing sin 2 6*jr % = constant ~ Q/c for a 
dipole geometry, we obtain 



pole 



(70) 



where fi m is the neutron star's magnetic dipole mo- 
ment. Substituting equation j70|) into (|69|), we obtain 



Jtot" 



(71) 



For the Crab pulsar, f^/z^/c 3 = 
fore, jtot ~ 1 is required, so that the observed 
spin-down luminosity 10 38 ' 65 ergs s" 1 may be realized. 
Analogous conclusions are derived for other rotation- 
powered pulsars. Moreover, the sharp pulse of the 
Crab pulsar may imply A po i <§; 7r(r* sin 6**) 2 ; there- 
fore, even j gap 3> 1 may be required. On these 
grounds, the limitation of jtot = jgap <C 1 derived 
for ji = J2 = were insufficient to apply to realistic 
pulsars. 

In the present paper, we relaxed the limitation of 
jtot by allowing ji or 32 to be non- vanishing. The 
results of the predicted 7-ray spectra are, therefore, 
more realistic compared with previous results ob- 
tained in Papers I, II, III, VI. However, even in this 
treatment, jtot is limited below unity. 

The next issue is, therefore, to consider whether 
we can construct an outer-gap model with super- 
Goldreich- Julian current density (i.e., j gap > 1). The 
Poisson equation (|30| ) tells that solutions exit even for 
ji + 32 3> 1, provided that j\ — 32 < 1- (For example, 
if ji = 32 1, the gap exists at the conventional 
null surface.) In this case, W becomes much smaller 
than those obtained for j\ + 32 < 1 because of the 
gap closure condition (eq. pof ). In the case of the 
Crab pulsar, the small W obtained for ji + j2 > 1 
fails the mono-enegetic approximation. To find solu- 
tions for jtot ~ ji + 32 ^ 1: we could assume much 
smaller collision angles so that the pair-production 
mean free path may become much larger. To settle 
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this issue, we must constrain the magnetic field ge- 
ometry around the gap and quantitatively infer the 
collision angles between the primary 7-rays and the 
secondary X-rays. 

In short, stationary gap solutions exit even for a 
super Goldreich- Julian current. In this case, the colli- 
sion angles should be much less than W/w-^c so that 
the emitted 7-ray flux may be consistent with ob- 
servations for the young pulsars whose X-ray field is 
dense (like Crab). For older pulsars whose X-ray field 
is less dense, on the other hand, we can in fact find 
solutions with super Goldreich- Julian current. It will 
be discussed in a subsequent paper. 

6.4. Comparison with Previous Works 

Let us compare the present methods and results 
with Paper V. In the present paper, E\\, N±(s), 
and G±(s,e 7 ) were solved from the Vlasov equa- 
tions for a non-vacuum gap, while in Paper V only 
-En field was solved from the Poisson equation for 
a vacuum gap, with the aid of W, which was de- 
duced from the gap closure condition. In the station- 
ary gap, the Vlasov equations automatically satisfy 
the closure condition; therefore, the obtained elec- 
trodynamic structures (e.g., W, Eu) are essentially 
the same between the two Papers, provided that the 
gap is nearly vacuum (i.e., jtot "C 1). By relaxing 
the boundary conditions of the magnetospheric cur- 
rent, and by solving the non- vacuum solution from the 
Vlasov equations, we first find in this paper an inter- 
esting behavior of the gap position: The gap shifts 
outwards (or inwards) when current is injected from 
the inner (or outer) boundary. The obtained GeV 
spectra are similar between the two papers, unless 
the gap is located well inside of the conventional null 
surface. In Paper V, the intrinsic TeV spectra were 
depicted in figure 6; on the other hand, in this pa- 
per, the TeV spectra after absorption were depicted 
in figures f§ |, |, |, and [l0| 

We briefly compare the present method with ZC97, 
who considered that the gap width is limited by the 
surface X-rays due to the bombardment of the parti- 
cles produced in the gap. The magnetospheric X-rays 
considered in this paper is much denser than the sur- 
face X-rays due to the bombardment. As a result, 
the localized gap in the present paper produces less 
intrinsic TeV flux compared with what would be ob- 
tained in ZC97 picture. 



6.5. Possibility of Another Solution Branch 

For cases 1, 2, and 3, the intrinsic TeV luminos- 
ity is comparable or less than the GeV one. There- 
fore, the Lorentz factors are limited primarily by the 
curvature process (eq.|y|). For case 4, however, the 
intrinsic TeV luminosity well exceeds the GeV one; 
therefore, the radiation-reaction forces are due to IC 
scatterings rather than the curvature process. In fact, 
we may expect a sufficient GeV flux via IC scatterings 
when the gap is located well inside of the conventional 
null surface. This is because the dense X-ray field will 
limit the particle Lorentz factors small (Paper II), and 
because the less-energetic particles scatter copious IR 
photons into lower 7-ray energies with large cross sec- 
tions (~ o"t)- There is room for further investigation 
on this issue. 
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